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SUMMARY 


A method is developed for solving the laminar and turbulent compressible 
boundary-layer equations for separating and reattaching flows. Results of 
this method are compared with experimental data for two laminar and three tur- 
bulent boundary- layer , shock-wave interactions. Several Navier-Stokes solu- 
tions were obtained for each of the laminar boundary- layer, shock-wave inter- 
actions considered. Comparison of these solutions indicates a first-order 
sensitivity in Cf to the computational mesh selected in both the viscous and 
Invlscid portions of the flow. ' 

Three turbulent, boundary-layer, shock-wave interactions were considered, 
one unsepafated interaction at M = 3 and two separated interactions at 
M = 1.47 and 3. Boundary-layer theory appeared to be adequate to describe the 
first two of these interacting flows. However, for the separated interaction 
at M = 3 , boundary-layer theory failed. 

Comparison of the present boundary- layer solutions with the Navier-Stokes 
solutions and with data for a given Mach number indicates, that as long as 
|vg/ug| is small, the boundary-layer approximation yields solutions whose 
accuracy is comparable to the Navier-Stokes solutions. A more general param- 
eter might indicate boundary-layer theory to be valid if some function of Mach 
number times lvg/ug| is small. 

Since the present boundary-layer solution is an Inverse method, that is, 
Cf(x) is specified rather than P(x),'a criterion for selecting the correct 
Cf distribution is required. It was originally anticipated that this cri- 
terion would be supplied by coupling the boundary-layer solutions to a shock- 
capturing, finite-difference solution to the inviscid-flow equations. 

Although several modes of coupling have been attempted, only slight success 
has been obtained to date because the shock-capturing scheme imposes a rela- 
tively thick shock wave on the viscous flow, which in turn requires an almost 
discontinuous change in flow angle at the boundary-layer edge. The inconsis- 
tency of this behavior in the two solutions has thus far obstructed attempts 
to obtain meaningful coupled solutions. Therefore, a truncated Navier-Stokes 
system of equations was examined and. it appears that this technique circum- 
vents these difficulties. 


INTRODUCTION 


It is generally agreed that the problem of separating and reattaching 
flows is one of the more challenging and technologically relevant problems in 
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computational fluid mechanics. In recent, years, the computational attack 
taken place on two fronts - the engineering approach by way of boundary- layer 
theory and the pure numerical approach by way of the Navier-Stokes equations. 
Reasonably complete summations of the status of these two approaches up^to 
about 1970 are contained in references 1 to 3. During the past two or three 
years, substantial additional progress has been made in the boundary-layer 
theory approach by Kllneberg and Steger (ref. 4) and Carter (ref. 5), and in 
the Navier-Stokes equation approach as represented by references 6 through 11. 
The basic position taken by each of these schools of thought can be summarized 
as follows: the Navier-Stokes equations are undoubtedly the correct equations 

for describing the flow fields in question and, since we can solve them, why 
settle for an approximation? The boundary-layer contingent would concur with 
the above but they would add that the cost of Navier-Stokes solutions both now 
and in the foreseeable future is far too high to permit their use as design 
tools. 

The present study examines these positions in some detail. In particular, 
we are concerned with establishing, at least qualitatively, the limits of 
applicability of the boundary-layer approximation for both laminar and turbu- 
lent separating and reattaching flows. Secondly, we are concerned with the 
development of an economical and reasonably accurate engineering calculation 
scheme for such flows. 

SYMBOLS 



C Chapman-Rubes in constant 

skin-friction coefficient, = t /(l/2)p u^ 
f ’ f w o o 

DC difference coefficient multiplying quantities at x 

DM difference coefficient multiplying quantities at x-Ax 

DP difference coefficient multiplying quantities at x+Ax 

E Eckert number, u^/2H 

F defined in equation (12) 

f dimensionless stream function 

H total enthalpy h + (u^/2) 

h static enthalpy 

a mixing length 

M initial Mach number 

o 
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number of nodes in x direction 

number of nodes in y direction 

pressure 

Prandtl number 

turbulent Prandtl number 

gas constant 

body radius 

entropy 

temperature 

defined in equation (12) 
velocity component in x direction 
velocity component in y direction 
streamwise variable 
cross stream variable 

weighting factor in type dependent differencing 

2C 

boundary-layer pressure gradient parameter = ~ 

ratio of specific heats 

difference in quantities 

boundary- layer thickness or stream angle 

displacement thickness 

residual error 


density 

defined in equation (18) 

viscosity or Mach angle (y = sln“^(l/M)) 


shear stress 



5 transformed x variable 

ri transformed ' y yariable 

( )' differentiation with respect to n 


• ' . ■■I 


Subscripts 


e 

evaluated at boundary-layer edge 

f 

evaluated at final 

conditions downstream of interaction 

i 

evaluated at ith 

X point 

j 

evaluated at jth 

y point 

£ 

lower boundary 


o 

. evaluated on initial line 

■t 

evaluated at stagnation conditions 

u ' 

upper boundary 


w 

evaluated at wall 

conditions 


ANALYSIS' ■ 


Viscous-Flow Calculations 

The procedure used to solve the viscous-flow portion of the calculation 
is a generalized Galerkin method (refs. 12 and 13) applied to the boundary- 
layer equations. The equations considered are 


9x 9y , 


3H . 3H 3 

^ ^ = 37 


3u . 3u dP . 3 / 3u . \ 

n O'' 37 ■ ■ 35? 37 "V 

lift » 

Pr 3y 2 \ Pr/ 


1\ 

3y 


:3y 


1 3H/3y . h 1 

Pr^ sS/iy ■ Pr^ 


P = pRT 


PT = 


3y 


3u 

3y 


( 1 ) 

( 2 ) 


(3) 

(4) 

(5) 
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Only adiabatic flows with Pr = Prj = 1 are treated here; hence equation (3) 
is replaced by h + (u^/2) = const. To suppress the streamwlse growth of the 
boundary layer and to facilitate initialization, equations (1) and (2) are com- 
bined under the Levy-Lees-Dorodnitsyn transformation (x,y) -*• (C»h) where 


2 k, PoV 


5 = f pupr dx, n = 
J ^o o o * 


m 


r ^dv 

Jo Po 




m- 


pu dy 


f = 


u 


” + If - H) ■ “ 


( 6 ) 


which results in 

(Cf) • + f f 

where Kj = t^2^/Pgr^ anH the primes indicate differentiation with respect to 
n« The following analysis closely parallels that of Kendall and Bartlett 
(ref. 13). The present method extends that of reference 13 to treat separated 
flow and is restricted to calorically and thermally perfect gases. 

y-dependent diffevenoing — To apply the generalized Galerkin method, the 
approximations for the stream functions, velocity, and shear between ad-jacent 
nodes are chosen as 


f = f + f’An + f" f f’’* Anf. + £**’ 

J+1 j + ^ 2 + 3 ^ ^j+1 24 

f * = f ' + f'An + f * ’ ' + f ' ’ ’ 

3+1 33^3 3 ^j+1 6 

fit _ fii . fill An^ ■ fill Ai]f_ 

3+1 ,33 2 ^ ^j+i 2 


(7) 

( 8 ) 
(9) 


These relations are obtained by a Taylor series expansion to terms including 

and substitution of f^^ = “ fj")/An. Equations (7) to (9) are 

substituted in equation (6) and the results Integrated with- respect to a unit 
square-wave weighting function from tIj to to yield: 


Cf 


^3+1 

- [2 + 2?E " 

’^j+l 

. + ff ’ 

+ KiPT 

"3 . 

"J -’fe' 

"3 


M+i 


"1 


'^3+1 9fi f^3+l 3f 

J, + 2C f" dn = 0 

^3 "^3 


(10) 
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where 


■* V 


and 


> f^lCPK = fIXPl + f''XP2 + fI*'XP3 + f!’’ XP4 = ff’^dn 

/ ^ j -j j j j+1 J 

k=i • 


XPl = An(fj 

> + f " + 

1 j 2 

\ 

XP2 = 


fi. 2An 

j 3 

XP3 = 


llAn 
j 15 



fu 4An 
j 5 


,, , , An_ ^ 

j 8 ^ ^ 
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The two terms remaining to be integrated are discussed under x-Dependent Dif- 
ferencing. This procedure, yields a consistent differencing method that is 
fourth-order accurate in n and equations (7) to (9) ensure that the stream 
function, velocity, and shear are defined ever3^here (not only at nodal points) 
and are continuous, everywhere. Equations (7) to (10) provide a sequence of 
4XN-1) equations in 4N unkno\ms, where- N is the number of nodes. normal. to 
the wall. Note that 3(N-1) of these equations are linear algebraic equations 
while (N-1) equations are nonlinear ordinary differential equations. The 
remaining four equations are Imposed as boundary conditions: 

i f = f ' = 0 at n = 0 


f = 1,. f", = 0 ■ at ,n = “ 

The complete system is solved by use of a Newton-Raphson Iteration scheme by 
differentiating with respect to . f ^ , f , . . . , f ’ " , f , and for sepa- 
rated flows 8» to yield 

'' j(f)Afj = ■ 

k 

where J(f) is the (4N x.4N) Jacobian, . e is the residual error, and Afj is. 
the incremental change in the solution parameters per Iteration. - 


To speed up the iteration process, the equations are ordered so that the 
3(N-1) linear algebraic equations plus the boundary condition f ’ (0) = 0 . 

occupy the first 3N-2 matrix rows. The matrix has the partitioned form (see 
ref. 14) : 


■ 

L2“ 




% 

_NLi 

Nl2. 




-^NL- 


The submatrix Li is a function of the nodal configuration only and can be 
inverted just once and used in all subsequent calculations. Formal manipula- 
tion and back substitution yield 

■ S‘[-=L - ^2“nl] 

The matrix ^NL 2 - Nl-iLj^L 2 j must be inverted every iteration, but it is only 
of order (N+2) . 

x-dependent differenoi.ng—'Ihe strearawlse differencing is carried out in a 
manner similar to that of reference 4. The term 



is decomposed and Integrated by parts as 


- f ^ + DCf^ + DMfJ^_Jdn = -[dp X fi+i jXPK + DC' f^ ^XPK 

J n . 


+DM £ f^_^ ^ jXPk] 


where, for- attached flow, 

DP = O’ , DC = 2/£n , DM = -DC ' ’ 

and for separated flow 

DP = (l-a)/£.n » a = 1 for f’ > 0,01 

DM = -(I4tx)/£n I a = lOOf , -0.01 f’ < 0.01 

DC = -(DP + DM) ; a = -1 , f ' < -0.01 


This procedure Incprporates backward differencing when the flow is in the. 
mainstream direction, forward differencing when the.; flow is reversed, and a 
central difference on and near the zero velocity streamline. The remaining 
term is decomposed as above. Integrated by parts, and expanded In Taylor 
series to yield; 


.. 157 



2? 


i: 


3f 

’j 


DC fif' 


'j+1 




'j+1 




• n. 


'j+1 


+ DC^f J jXPK^ 


+DPi:f^ + Dm2>J XPK^_^ 


where, for attached flow, excluding the points of separation and reattachment 
yields 

DC = 2/ln , DP = 0 , DM = -DC 

and for separated flow, 

DC = 0 , DP = 2/ln , DM = -DP 


As shown, a straightforward marching routine Is used In the attached flow 
while, for separated flow, the entire separation region must be relaxed 
simultaneously. 

The present method requires no under relaxation and a typical \»ell- 
separated flow with 20 n points and 10 X points can be relaxed to a maximum 
residual 0(10“**) In about 22 Iterations compared to 800 for the method of 
reference A and 100 to 200 Iterations for the method of reference 5. Although 
the present method may require more operations per iteration, it is at least 
competitive with other methods cited. The speed of the present method 
primarily results from the spllned Taylor series used as approximating func- 
tions. First, they provide an effective fourth-order-accurate, finite- 
difference representation in the y direction which permits a relatively 
sparse nodal array to yield high accuracy and, second, the spline character 
enhances the stability markedly. 


Inviscid-Flow Equations and Coupling Schemes ... 

In this section, the technique used to solve for the Invlscld flow field 
and to obtain the effects due to the boundary layer on that flow Is described. 
A technique for obtaining fully coupled interactive solutions of the Invlscld 
and ylscoiis flows is being developed, but is not sufficiently advanced to ' 
discuss here. Rather, at this point, boundary-layer solutions for a pressure 
and skin-friction distribution, appropriate to the data of reference 13, have 
been obtained and then various schemes for, matching this solution with the • 
invlscld flow have been investigated. The invlscld flow is for supersonic 
flow between a shock-wave generator and some matching line given by the 
boundary-layer solution. 

Inoisaid-flow equations — The basic equation used In the Invlscld flow- 
field analysis Is similar to that of reference 15. The key features, however. 
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are described here for completeness. ■ In a Cartesian coordinate system, the’ 
steady, inviscid, two-dimensional fluid' dynamic equations (continuity, x and y 
momentum) are written in conservative form as 

U + F = 0 (11) 

X . y 

The three components of the vectors U and F, which represent the conservative 
variables, are defined as 



pu 


'pv 

u = 

kp+pu^ 


puv 


puv 


kp+pv^ 


( 12 ) 


where k = (y~1)/2y. The units of these equations were normalized by dividing 
both pressure and density by their respective stagnation conditions, while the 
velocities were divided by the maximum adiabatic velocity. With this normaliz- 
ation and the further restriction that the flow be adiabatic, the energy equa- 
tion can be written: 

^ + u2 + v2' = 1 (13) 

. p ) 

To calculate two-dimensional flow between two nonparallel walls, in this 
case the shock-wave generator and the lower coupling boundary, it is most con- 
venient to normalize the coordinate system so that the upper and lower bound- 
aries become parallel. The transformed coordinates are 



where y^:= Yj, = Y£,(x). Applying this transformation to. equation 

(11) results in 


U + F^ = 0 (15) 

X C 

where U = CU arid F = F + u[y£(5-l) - yu"| • this section, the primes denote 
differentiation with respect to x.) 

The above differential equation is, integrated using MacCormack's second- 
order-accurate differencing scheme (see ref. 15). For flow-field points, that 
is, points on neither the upper nor lower boundary, the predictor and corrector 
equations are _ , 

. “r‘ ' - (proctor) ( 16 ) 

Subscript j identifies particular points on a data line as follows: 1, 
lower body point; 2, 3, . . . , flow-field points; and N^, upper body 

point. . ,The data lines are Identified by the superscripts n and n+1. The 
bars indicate a predicted value. 
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(17) 


f' 


% 




(corrector) . 


The above equations are applied to the three vector components of U. 




The step size in the x direction is found by determining the maximum 
slope of the characteristic surface, inclined at the local Mach angle p rela- 
tive to the local stream angle 6, for each point as follows: 


Oj(±) = [tan(±w + 6)"| 


(18) 


where the ± denote positive and negative directions of the Mach angle from 
each point. The entire data line n is surveyed to determine |g^(t) ' 

Depending on the sign associated with the maximum slope, Ax is found from 


or 


Ax(+) 


Ax(-) 


pj+1 

- 


n+r 

L 

f”n n 

_n+l n+1 

n+1 


h' 

- 



L 

n. . 

max 




(19) 


( 20 ) 


where K is a constant that controls the step size. If K = 1, the CFL con- 
dition is satisfied identically. 


Application of the finite-difference equations (16) and (17) yields pre- 
dicted and corrected values of the three components of the conservative- 
variable vector U on the new data line. These must be decoded after the 
predictor step and corrector step in turn. This decoding into physical vari- 
ables is accomplished by solving the following: 


V 


u 


P 



n+1 

i 


U 2 + - 4U2k(l-k)(l-v2) 

— 2Ui (1-k) 


Cui 


n+1 


P = pCl-u^-v^) 


n+1 

j 


( 21 ) 


/ 


Constructing a technique that will give a physically realistic solution 
of the conditions on a solid boundary is difficult in a predictor-corrector 
sequence since strict application requires an imaginary point inside the body. 
Abbett (ref. 16) developed a technique that satisfies the surface tangency 
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condition and relies on information provided by the finite-difference equa- 
tions. For the upper and lower surfaces, respectively, the predictor equations 
are written as 


and 





= uf - || (F2 - Fi) 


( 22 ) 

(23) 


Modified corrector equations that maintain second-order accuracy at the 
upper and lower surfaces, respectively, are given by 

C' = ik + - 1 - 2 >] <“) 


and 


uf ' . + uf * - ^ (rf ‘ - Ff ') - H (F? - 2F? - F?)l 


(25) 


After the corrector equations are decoded at the boundaries, the flow 
variables will not necessarily satisfy the boundary conditions at the wall. 
The technlq-ue for satisfying the boundary conditions depends on a scheme for 
matching the boundary-layer solution. 

The coordinates of the upper boundary are determined from 


y = y + X tan 6 (26) 

•^u ■^u u 

o 

where 6^ is the deflection angle of the upper wall. When matching to the 
output of the boundary-layer solution, the lower boundary is specified by a 
table of coordinates. These coordinates are used to fit a parabola of the 
form 


y = ax^ + bx + c (27) 

between each set of three consecutlye points. The three points used to deter- 
mine the coefficients in the above equation were chosen so that two of the 
points should have x > x”"*"! . The slope of the boundary is found from the 
derivative of equation (27) . 


Matching scheme — Two techniques for matching the inviscld and viscous 
solutions were investigated: (1) requiring tangency along the matching con- 

tour and (2) forcing the inviscid solution to agree with the u and v veloc- 
ity components from the viscous solutions along the matching line. For the 
first matching scheme, the flow angle from the decoded corrector equations on 
the boundary is foimd from 


6 = tan~^ — (28) 

u 
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This flow angle is compared with the slope of the boundary at . If the 

two angles disagree, then a Prandtl-Meyer turning is applied at the new body 
point using 

* W . 

= 1 - Y^A6 (M^) - i]a62 (29) 

where A6 is the difference in the two angles and 

Af = (30) 

Note that the sign of A6 will result in either an expansion (A6 > 0) or an 
isentropic compression (A6 < 0) and the pressure, resulting from equation (19), 
is taken as the predicted pressure on the body at data line N+1, The density 
is found from 

= [pcp^'' (31) 

-s-s^/V . . 

where C = e and S^. is some reference entropy. Once p and p are 

known at the boundaries, the total velocity is obtained from 

q = Vl (p/p) ' ' (32) 

The velocity components then result from the tangency condition at the. bound-: 
aries. This technique was applied on the upper boundary throughout both 
matching schemes investigated. With this matching scheme, two different 
matching boundaries were considered - 6* and an arbitrary streamline where the 
Mach number remained supersonic throughout 'the interaction. Results are pre- 
sented only for the latter case since similar results were obtained for both 
matching contours. 

The second matching scheme, forced the inviscld solution to match values 
of u and v given by the viscous solution along the matching contour, here 
taken as the boundary-layer edge. This technique does not Impose the tangency 
condition along the matching boundary, but rather allows for mass exchange 
from the inviscid to viscous portions of the flow, a physically realistic 
situation. Since u and v along the, matching contour are now taken to be , 
known, these values are used in the decode; equations (21) to solve for ' p^+l 
anci’ . , . ■ 


Navier-Stokes Solutions 

Solutions to the Navier-Stokes equations presented here were obtained 
from several sources. Solutions of MacCormack (ref. 3), MacCormack and Bald\d.n 
(ref. 11), Mess^lna (unpublished), and Skoglund and Gay (ref. 6) were obtained 
from the sources cited, while additional solutions were obtained in the. pres- 
ent study using the computer codes of MacCormack and Baldwin (ref. 11) and 
Carter (ref. 10). 


162 


Space does not permit a detailed description of the differences and simi- 
larities of the several codes considered. There are, hov/ever, some striking 
differences and similarities in the several solutions presented x^hich warrant 
discussion in some detail. Figure 1 shows the envelope of four Navier-Stokes 
solutions obtained from references 3 and 10 to 12, together with the data of 
Hakkinen et al. (ref. 17), for an unseparated, laminar-boundary-laver , shock- 
wave interaction at M = 2. Two points immediately come to mind: first, the 

large discrepancy in skin friction and the relatively small discrepancy in sur- 
face pressure. At first glance, this implies a low sensitivity of the pressure 
distribution to the skin-friction distribution and a lack of uniqueness in the 
Navier-Stokes solutions. However, 'a sequence of solutions obtained here using 
MacCormack and Baldwin's (ref. 11) code, with successive mesh refinement in 
the outer inviscid flow, shows differences in Cf of the same order as those 
shown in figure 1. This arises from the fact that for the coarsest mesh con- 
sidered, the externally generated shock wave takes on a thickness that, x^hen 
projected onto the boundary layer, is of the same order as the length of inter- 
action. The viscous flow is then responding to a continuous compression 
rather than an imposed shock wave. As the nodal structure was refined, the 
skin-friction distribution approached a single curve. 

Figure 2 shows a similar envelope of Navier-Stokes solutions obtained 
from references 6, 8, 10, and 11, together xjith experimental data for a well- 
separated, laminar-boundary- layer , shock-wave Interaction (ref. 17). The 
results here are similar to those discussed above and the same conclusions can 
be draxm. 

The point of the foregoing presentation is not to show how poorly Navier- 
Stokes solutions perform but rather to induce the user of Navier-Stokes codes 
to examine his results critically and, in particular, to examine the computa- 
tional mesh dimensions in the light of the smallest relevant physical scale in 
the problem under consideration. 

In subsequent sections, the boundary-layer calculation scheme developed 
is compared to the Navier-Stokes solutions of MacCormack (ref, 3). 

• . . I 

RESULTS 


The results of the inverse boundary-layer method are compared with 
Navier-Stokes solutions and xd.th experimental data. The flows considered all 
occur on a flat plate and the pressure rise is caused by an externally gen- 
erated shock wave. The specific flow parameters and experiments considered 
are listed in table 1. 


Table 1 


Reference 

M 

00 

Re 

Pf/Po 

Remarks 

17* 

2 

1.84x10=’ , 

1.2 

Laminar 

17* 

2 

1.96x105 

1.4 

Laminar 

18 

1.47 

4xl05/ft 

2.25 

Turbulent 

19 

2.93 

5.7xl0'^/m . 

2.5 

Turbulent 

19* 

2.93 

5.7xl0^/m 

5.0 

Turbulent 


*Navier-Stokes solutions available. 
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Finally, we will describe the progress made, to date, in coupling the inverse 
boundary-layer method to the inviscid calculation scheme to provide a complete 
flow-field prediction method. 

. Laminar-Boundary-Layer, Shock-Wave Interaction 

Two laminar-boundary-layer, shock-wave interactions are considered; an 
unseparated interaction at M = 2.0 having an overall pressure rise Pf/Po . 

1.2 and a well-separated interaction at M = 2.0 having an overall pressure 
rise Pf/Po 1.4. These two interactions correspond to the experiments of 
Hakkinen et al. (ref. 17), as shown in figures 6a and 6b, respectively, of 
that report. For convenience, these interactions will be called Hakkinen 6a. 
and 6b. 

Note that in this and all subsequent comparisons v;ith Navier-Stokes solu- 
tions, the skin-friction distribution obtained from the Navier-Stokes solution 
was input to the inverse .boundary-layer method. This procedure ensures that a 
comparison of other parameters, that is, velocity and pressure distributions, 
provides a true measure of the validity of the boundary-layer assumption, 
without the peripheral considerations of downstream boundary conditions or 
matching conditions between viscous and inviscid flow. However, some differ- 
ences can arise even in this case because of the order of the difference 
approximation used. The present inverse method is third-order accurate in 
wall shear as opposed to first-order accuracy imposed by assuming a linear 
variation of velocity with distance away from the wall in the Navier-Stokes 
solutions. An examination pf the differences attributable to this assumption 
indicates:- that Cf . can vary ±10 percent for the same velocity profile, depend- 
ing, on the order of curve fit used to deduce the derivative at the wall. 

Figure 3 compares the pressure distribution obtained from present 
methods with a refined mesh Navier-Stokes calculation using Hung and 
MacCprmack’s code (ref. 9) for Hakkinen 6a. The two solutions are nearly 
indistinguishable everywhere and they. are in reasonable agreement with the 
data. In figure 4, the velocity distributions at the x location correspond- 
ing to the minimum value of Cf (fig. 3) are compared. The differences 
between the present method and MacCormack’s solution are not large and are 
largely attributable to differences in transport properties, that is, Pr = 1 
and y/Uo “ '^/'^o present method and Pr = 0.72 and y from the 

Sutherland’ law- in MacCormack’s code. 

Figure 5 shows the streamwlse pressure distribution for- Hakkinen 6a, for 
y = 0, 6/2, 6, as computed by MacCormack (ref. 3). These curves are essen- 
tially indistinguishable and hence, within the boundary layer, 3p/3y z 0. 

This, in turn, implies a trivial solution to the y momentum equation which, 
from the usual boundary order-of-magnitude arguments, indicates that the 
boundary-layer equations are adequate for this flow. 

Figure 6 compares the pressure distribution obtained from the present 
method with that of MacCormack' s solution for Hakkinen 6b. For this case, a 
slight difference in surface pressure distribution is obtained near x/L = 0.15. 
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Figure 7 compares the velocity profiles obtained from the present method and . . 
from MacCormack's method for Hakkinen 6b. Differences between the profiles, 
are explicable in terms of differences in transport properties discussed 
previously. Figure 8 shows the strearawise pressure distributions for Hakkinen 
6b, as computed by MacCormack (ref. 3). This figure shows the normal pressure 
gradient developing in about the same region noted in figure 6 (x/L =0.15 
corresponds to x/Ax = 20). This normal pressure gradient appears to be an 
Inertial effect associated with strong streamline curvature in the outer por- 
tion of the boundary layer where M(y) -»■ Mg. For this case, the boundary- 
layer method predicts a relatively large value of [vg/Ugl (not shown), which 
indicates the beginning of a breakdown of the boundary-layer approximation. 


Turbulent Boundary-Layer, Shock-Wave Interaction 

The question of the validity of boundary-layer theory in turbulent 
boundary-layer,' shock-wave interactions is somewhat less clear than for lami- 
nar interactions. This is true for three reasons. First, few Navier-Stokes 
solutions with simple turbulence models are available for comparison (ref. 18). 
Second, the current turbulence models are demonstrably Inadequate for flows 
with rapidly changing boundary conditions and, third, only a few experiments * 
are available in supersonic turbulent flow which are sufficiently detailed to 
permit a useful comparison. Because of these facts, only a cursory comparison 
of results of the present method with both Navier-Stokes solutions and experi- 
mental data can be made. 

Figure 9 compares the present method with the data of Seddon (ref. 19), 
for a normal shock-wave, turbulent-boundary-layer interaction with an initial 
Mach number of 1.47. For this case, the experimentally determined value of 
Cf(x) was input to the inverse, boundary-layer method in the attached-flow 
regime. For the reversed-f low region, where no data were available, three 
estimates of were made (defined by the solid, dashed, and dot-dashed 

lines in the upper half of figure 9) . The corresponding pressure distributions 
are shown in the lower half of figure 9. Two conclusions can be drawn from 
this figure: first, the boundary-layer equations are capable of reproducing 

the observed behavior and, second, the pressure distribution is sufficiently 
sensitive to variations in Cf to permit a coupled solution to distinguish 
between similar distributions of skin-friction coefficient. 

Figure 10 shows a comparison of a direct calculation using the present 
method. with the data of reference 20. The experimental skin-friction distribu- 
tion shown was deduced from the measured mean velocity profiles using the 
method of reference 21. These data were taken at a free-stream Mach number of 
2.93 and a unit Reynolds number of 5.7x10^ m”^. The pressure rise was effected 
by a shock generator set at 7° incidence to the oncoming flow. As noted above, 
the predicted results were obtained for this case by calculation in the direct 
mode, that is, P(x) was specified and Cf computed. The two predicted skin- 
friction distributions shown in the upper portion of figure 10 result from two 
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turbulence models used for this calculation. The equilibrium model is that 
described above while the exponential lag model is that suggested by Rose and’’" 
Johnson (ref. 22) and subsequently used with some success in the Navier-Stokes ^ 
calculations of reference 23 (see also ref., 18). One may conclude then that,*'" 
given the appropriate turbulence model, the boundary-layer equations are ade-^ 
quate to describe flows of this type. 

The last case considered was also taken from references 20 and 21. The 
initial flow properties were the same as those discussed in figure 10 but the 
shock generator angle was set at 13° to the oncoming flow. As shovm in the 
upper portion of figure 11, the flow has an extensive region of separation 
(x - X ) ~ 46,. For this case, the inverse boundary-layer method ; 

failed. The predicted pressure distribution shown in figure 11 Indicated . - 
nearly twice the pressure rise observed experimentally. Direct solutions 
employing the experimentally observed pressure distribution failed to predict * 
any separation. This might be attributed to three-dimensional effects in the 
experiment except that the Navier-Stokes solutions of Baldwin and Rose (see 
ref. 18), using a similar turbulence model, did a plausible job of predicting 
the flow. An examination of this solution shows extensive regions of signif- 
icant normal pressure gradient. 


Coupled Solutions 


Figure 12 compares the pressure distributions resulting from the two 
coupling schemes for the Hakkinen 6b case with that given by the Navier-;- ^ . 

Stokes program of reference 3. Neither coupling scheme appears to' be ade- 
quate when applied in a noniterative manner. The tendency of shock-capturing 
techniques to smear the shock wave over several grid points and to overshoot 
the pressure rise through the shock-wave reflection is evident. Perhaps the 
best overall qualitative agreement with the Navier-Stokes results is provided 
by. the u-v coupling scheme; however, as noted earlier, the discontinuity in 
V that results from the boundary-layer calculations will strongly affect this 
scheme. 
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Truncated Navier-Stokes Equations 


To circumvent the difficulties encountered in coupling viscous and 
inviscid flows, an alternative procedure was considered. Rather than solve 
coupled viscous and inviscid equations in an iterative mode, we proposed to 
solve the following system of equations: 
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together with an equation of state and a consistent energy conservation rela- 
tion. This was accomplished by simply removing the appropriate terms from the ' 
computer code of MacCormack and Baldwin (ref. 11). These equations converged 
everywhere within 1 percent pf full Navier-Stokes solution for the Hakkinen 
cases 6a and 6b using the same mesh configuration (results not shown). The 
above method was slightly more efficient than the solution of the full Navier- 
Stokes equations because fewer operations need to be performed in the computer 
code. 

One can deduce, however, that these equations can be solved and since 
they are parabolic in the streamwise variable, they can be solved by ai forward- 
marching procedure. This latter fact would allow a first-order improvement in 
computational efficiency. Unfortunately, these conclusions were reached very 
late in the present study and time did not permit significant exploitation of, . 
the properties of these equations. 


CONCLUDING REMARKS 


We have demonstrated the validity of the boundary-layer equations applied 
to certain shock-wave, boundary-layer interactions. As has been suggested 
(e.g. , ref. 2), solutions to the boundary-layer equations begin to depart sig- 
nificantly from those of the Navier-Stokes equations under the same conditions 
for which nontrivial normal pressure gradients are first observed in the 
Navier-Stokes solutions. Based on the physical argument that high-velocity 
flow, turning through a large angle, requires the action of a large force over 
a short distance, ■ one may deduce that the normal pressure gradient arises pri- 
marily from the Inviscid characteristics of the flow field. From this one may 
conclude that when some parameter involving both Mach number and turning angle, 
for example, Mg|vg/Ug|, exceeds a critical value, boundary-layer theory will 
fail. 

Despite the fact that boundary-layer theory can be shown to fail at high 
supersonic Mach numbers when extensive separations occur, it would appear 
that, for many technological problems, boundary-layer theory is quite satis- 
factory. This should be particularly true in the design of transonic airfoils 
and for engine Inlets to be used on transonic and low supersonic flight 
vehicles when high angle-of-attack trajectories are not required. 

Note that a typical solution presented here required about 3 min of CPU 
time on the IBM 360/67. Our experience indicates a factor of about 25 reduc- 
tion in CPU time when the CDC 7600 is used. While the present Inverse method 
is not yet a complete predictive method, it shows substantial promise if the 
coupling problems can be solved. Alternatively, if these problems defy solu- 
tion, the truncated Navier-Stokes system described earlier appears to be 
potentially useful. 
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Figure i«- Comparison of envelope of four, Navier-Stokes solutions with the 


data of Hakklnen et al". ; H » 2 , Re 


shock 


2.84x105, Pf/pQ = 1.20. 



Figure 2.- Comparison of, envelope of four Navier-Stokes solutions with the 
- data of Hakklnen. et al.; M * 2, Re . v , = 2.96xio5 p^/p = ,1.40. 
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Figure 5.- Predicted pressure distributions for an unseparated laminar 
boundary-layer, shock-wave Interaction; calculations of JlacCormack. 
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Figure 6.- Comparison of results of present method with calculations of 

MacCormack. 
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Figure 7.- Comparison of predicted velocity profiles at minimum C^. 
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Figure 8,- Predicted pressure .distributions for a separated laminar Jbioimdary 
layer, shock-wave interaction; calculations of MacCormack. 
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Figure 9.- Comparison of fhe results of the present method with the data of 

Seddon (ref. 18). 
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Figure 10.- Comparison of the present method with the data of reference 19; 

turbulent unseparated flow. 
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Figure 11.- Comparison of the present method with the data of reference 20 

turbulent separated flow. 



Figure 12.- Comparison of coupling schemes. 


